The scattering problem in non-equilibrium quasiclassical theory of metals and 
superconductors: general boundary conditions and applications 
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I derive a general set of boundary conditions for quasiclassical transport theory of metals and 
superconductors that is valid for equilibrium and non-equilibrium situations and includes multi-band 
systems, weakly and strongly spin-polarized systems, and disordered systems. The formulation is 
in terms of the normal state scattering matrix. Various special cases for boundary conditions are 
known in the literature, that are however limited to either equilibrium situations or single band 
systems. The present formulation unifies and extends all these results. In this paper I will present 
the general theory in terms of coherence functions and distribution functions and demonstrate its 
use by applying it to the problem of spin-active interfaces in superconducting devices and the case 
of superconductor/half-metal interface scattering. 

PACS numbers: 74.20.-z, 74.45. +c, 74.81.-g 



I. INTRODUCTION 

For the theoretical understanding of transport in met- 
als and superconductors Landau's concept of quasipar- 
ticles acting as elementary excitations over the ground 
state has been of immeasurable valuei^ In a normal 
metal, electrons are in a strongly quantum correlated 
state due to Pauli's exclusion principle and due to 
Coulomb interactions. Conduction electrons in metals 
are, however, quasiparticles, i.e. elementary excitations 
in the vicinity of the Fermi surface that are rarely scat- 
tering with each other as a result of phase space restric- 
tions. Although these quasiparticles are strongly coupled 
to electrons far away from the Fermi surface, renormal- 
izations due to these interactions are constant over the 
energy range of interest {k^T, with temperature T) and 
thus can be treated as phenomenological parameters of 
the theory* 1 * 2 - Quasiparticles are represented by a classical 
distribution function and obey a semiclassical Landau- 
Boltzmann transport equation. 1 

Landau's Fermi liquid theory can be formulated in 
a systematic way within a diagrammatic expansion of 
many-body Green's functions* 3 - Asymptotic expansion in 
the small parameter k-QT/Ep (with the Fermi energy Ep) 
leads to the quasiclassical theory of metalsj 4 * 5 -^ that de- 
scribes the range (kp,T) 2 /Ep <C kp,T <C E-p in tem- 
perature well. In leading order, the dynamical equa- 
tions for Green's functions can be transformed into Lan- 
dau's transport equation for quasiparticle distribution 
functions Electrons that are far away from the 
Fermi surface and thus do no represent quasiparticles en- 
ter this theory as effective interaction vertices. Only a 
small number of these vertices is needed to describe the 
dynamics of the quasiparticles. 

The development of semiclassical concepts for the su- 
perconducting state was pioneered by Geilikman 8 -^ and 
Bardeen et alr^ soon after the development of the BCS- 
theory of superconductivity,!! Several early work o 12 ' 13 : 14 
on transport and linear response in superconductors 



showed that various semiclassical concepts of Landau's 
Fermi liquid theory could be readily generalized to the 
superconducting state. A formulation of the equilibrium 
theory of superconductivity near the superconducting 
critical temperature T c in terms of classical correlation 
functions was developed by de Gennes,! 5 - 

In the seminal works of Larkin and Ovchinniko\*!£ and 
Eilenberger— the concepts of the BCS pairing theory 
of superconductors 1 - 1 - were merged with the concepts of 
Boltzmann transport equations within Landau's Fermi 
liquid theory. This quasiclassical theory of superconduc- 
tivity was later generalized to non-equilibrium phenom- 
ena by Eliashberg^ and Larkin and Ovchinnikov. 18 

Quasiclassical methods can be applied to both wave- 
function techniques and Green's function techniques. In 
the former case, the starting point are Bogoljubov's 
equations ; 15 : 19 leading in quasiclassical approximation to 
Andreev's equations for the envelopes of the waves, 2 ^ Al- 
ternatively, one can start from the microscopic Nambu- 
Gor'kov matrix Green's functions, 21 - In quasiclassical 
approximation they result into envelope Green's func- 
tions that vary on the coherence length scale, £o — 
h\vp\/2irkp,T c (with Fermi velocity vp), and the time 
scale to = h/A (with gap A), and are free of irrelevant 
fine-scale structures on the Fermi wave length scale. 

Dynamical phenomena are described within quasi- 
classical theory by using the Keldysh Green's function 
technique. 22 Quasiparticle states in superconductors are 
coherent mixtures of particle and hole states. The degree 
of mixing is determined by the superconducting order pa- 
rameter A. The spectrum of quasiparticles is coupled to 
quasiparticle distribution functions, and this coupling is 
expressed in Keldysh's technique by two types of Green's 
functions, g R,A and g K , that are elements of a 2x2 ma- 
trix g. The information about distribution functions is in 
the Keldysh part, g K . Different formulations in terms of 
dynamical distribution functions in the superconducting 
state have been introduced by Larkin and Ovchinnikov, 18 
by Shelankovj^ 3 - and by the author^ 
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The derivation of boundary conditions for quasiclas- 
sical Green's functions is a difficult problem. For mi- 
croscopic Green's functions the formulation of bound- 
ary conditions, e.g. in terms of scattering matrices or 
transfer matrices at interfaces, is rather simple. In con- 
trast, in quasiclassical theory only the envelope function 
of the Bloch waves is known. The information about the 
phase of the waves is, however, missing. Under these 
circumstances it is not a priori clear if boundary con- 
ditions can be formulated within quasiclassical theory. 
That this is indeed the case was shown independently 
by Shelankov^ and by Zaitsevi 2 ^ More general formula- 
tions have been derived subsequently ; 27 ' 28 ' 29 ' 30 including 
a formulation in terms of scattering matrices by Millis, 
Rainer and Sauls. 30 However, owing to the normalization 
condition for the quasiclassical propagator, the boundary 
conditions so far were formulated as non- linear equations. 
Furthermore, their practical use was limited as they con- 
tained unphysical, spurious solutions that lead to insta- 
bilities in numerical calculations. 

Progress has been achieved by using the projector for- 
malism of Shclankov, 23 that allows an explicit formula- 
tion of boundary conditions for both equilibriu m 31 ' 32 ' 33 
and non-equilibrium 3 - 2 - situations. These boundary con- 
ditions have been generalized for the single band case 
to include spin-active interfaces in equilibrium 3 ^ and in 
non-equilibrium, 35 diffusive interface scattering, 36 and 
interfaces with strongly spin polarized ferromagnetsi 37 ' 38 
An alternative, equivalent, route has been followed via 
transfer matrices! 39 ' 40 ' 41 ' 42 ' 43 All the developments above 
were complemented by boundary conditions for diffusive 
superconductor a 44 ' 45 ' 46 that are appropriate for the dif- 
fusive limit of quasiclassical theory, the Usadel theory^ 7 . 

In this work, we will pursue the approach in terms of 
scattering matrices, and will present the boundary condi- 
tions in their most general form. Our results include all 
previous formulations as special cases, and are capable of 
describing e.g. non-equilibrium effects, multiband met- 
als, spin polarized systems, and diffusive interfaces. In 
most of these cases the present formulation leads to more 
transparent and compact boundary conditions, that al- 
low (i) for a very effective numerical implementation and 
(ii) better analytical treatment due to their simpler struc- 
ture. We use throughout the notation of Ref. [32|. 



II. THEORETICAL DESCRIPTION 

Quasiclassical theory is a powerful tool for describ- 
ing inhomogeneous superconducting systems in and out 
of equilibrium, covering both ballistic and diffusive 
materials^ii 7 .^^^^^^ 2 . All the relevant physical 
information is contained in the quasiclassical Green func- 
tion g(e, p F , R, t). Here e is the quasiparticle energy mea- 
sured from the chemical potential, p F the quasiparticle 
momentum on the Fermi surface (that can have several 
branches), R is the spatial coordinate, and t is the time. 
The "hat" refers to the 2x2 matrix structure of the prop- 



agator in the Nambu-Gor'kov particle-hole space, and the 
"check" to the 2x2 Keldysh matrix structure. The equa- 
tion of motion for g is the Eilenberger equation ; 16 ' 17 

[efsi — h,g] + ihv F -'Vg = (1) 

subject to the normalization condition 

gog = -7T 2 1. (2) 

The elements of the 2x2 Keldysh matrices are matrices 
in Nambu-Gor'kov particle-hole space, 

»-(££)■ *-(o"E)- (3) 

and T3 is the third Pauli matrix in particle-hole space. 
The o-product combines a time convolution and a ma- 
trix product and is explained in Appendix [AJ For what 
follows it is useful to think about it as discretized in time, 
in which case its properties are that of conventional ma- 
trix multiplication^ 3 - In equilibrium we will have to re- 
tain a matrix structure if the spin degree of freedom is 
active, in which case the o-product reduces to a matrix 
multiplication in Pauli spin space. 

Self energies enter Eq. (fT]) via the matrices 

where diagonal (S) and off-diagonal (A) self energies are 
determined by self-consistency equations. In this paper 
we do not, however, need to specify the exact form of 
these equations, and will assume for what follows that 
their solutions are given. There are fundamental symme- 
tries that relate the particle and hole components of both 
self energies and Green's functions £ We express these 
symmetries throughout this paper by using the particle- 
hole conjugation operation that is defined in the mixed 
(e, t) representation via 

Q(z,p F ,R,t) = Q(-z*,-p F ,R,t)*, (5) 

where z — e is real for the Keldysh components, and z is 
situated in the upper (lower) complex energy half plane 
for retarded (advanced) quantities. 

The characteristic curves of the partial differential 
equation (TTJ) define the quasiclassical trajectories. Tra- 
jectories are labeled by the position on the Fermi sur- 
face, p F , and are aligned with the Fermi velocity v F (p F ). 
Quasiparticles move along these trajectories, thereby be- 
ing coherently coupled to the condensate. 

Eqs. ([T]) and @ must be supplemented by boundary 
conditions at the two ends of each trajectory. Eq. (TTJ) is 
numerically stiff, with exponentially growing solutions in 
both directions. In addition, unphysical solutions must 
be eliminated using the normalization condition Eq. 
Both problems are solved in a natural way with the pa- 
rameterization of the quasiclassical Green's functions by 
coherence and distribution functions.- 3 ^ These are phys- 
ical quantities that obey initial value problems with a 
stable integration direction, and automatically ensure the 
normalization of g. 
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FIG. 1: (Color online) a) The coherence function y(t,t') de- 
scribes the local probability amplitude for conversion of a hole 
(dotted line) at time t' to a particle (full line) at time t. For re- 
tarded functions t > t' , and for advanced functions t < t' . b) 
The corresponding local amplitude for conversion of a particle 
at time t' into a hole at time t is described by the coherence 
function 7(4, t'). c) Diagrammatic representation of Eqs. l[7| 
and (fTTj) . 

A. Coherence functions and distribution functions 

The numerical solution of the (non-linear) system of 
Eqs. (P) and @ is greatly simplified by using a conve- 
nient parameterization of the Green functions in terms 
of retarded and advanced coherence functions 7 R A , 7 R A , 
and distribution functions x, £ i 24 i 32 i 54 i 55 The coherence 
functions are a generalization of the so-called Riccati 
amplitude o 56 ' 57 to non-equilibrium situations. Using a 
projector formalism as described in Appendices [B] and 
[C]we can write the retarded and advanced Green's func- 
tions [here the upper (lower) sign corresponds to retarded 
(advanced)] as 

9 R A = T2m (_ Q ~ _JT ±wrf 3 , (6) 

with the parameterization 24 

Q = (l-7°7)~\ F = {1 - l°iy lo l, (?) 

g = (1-707)-!, f= (1 - 707)^07. (8) 

The inverse is denned via the o-product, 

(...)- 1 o(...) = (...)o(... )-i = l, (9) 

with the unit element 1 (see Appendix [X| . Obviously, 
we can calculate the coherence functions from 

1 = g- 1 oT = Tog-\ ^ 1 = g- 1 oT = Tog- 1 . (10) 

In order to obtain a diagrammatic representation we re- 
formulate the problem in terms of Dyson equations 

g = 1+^/0707, ^ = 7 + ^0707, (11) 
g = 1 + £0707, T = 7 + ^0707. (12) 



In Fig. [T] the corresponding diagrammatic expansion is 
shown. Here, and in the following, we adopt and ex- 
tend a diagrammatic notation by Lofwander, Zhao and 
Saulsj 58 ! 59 ' 60 The quantity g describes the local spectral 
amplitude of a particle-like excitation in the presence of 
a condensate. This amplitude is renormalized from its 
normal state value S(t — t') due to multiple virtual An- 
dreev scattering processes that take place in the presence 
of an off-diagonal complex condensate field A. The same 
holds for hole-like excitations, described by the quantity 
g. The "anomalous" propagators T and T result from 
the local coherence amplitudes for particle-hole conver- 
sion, 7, and for hole-particle conversion, 7, again by tak- 
ing into account renormalization due to multiple virtual 
Andreev processes. For small superconducting ampli- 
tudes (e.g. near T c ) the anomalous propagators coincide 
with the coherence amplitudes. The four functions T, 
g and g are inter-related via 7 and 7, and a number of 
identities hold that are shown in Fig.[2]diagrammatically. 

Although the coherence functions 7 and 7 are suf- 
ficient to describe the retarded and advanced Green's 
functions, the quantities in Eqs. pip - (fT2]) allow for an 
effective formulation of boundary conditions (see below). 
The Keldysh part of the propagator can be formulated 
in terms of these and a suitable distribution function for 
particle-like and hole-like excitations, respectively, in the 
following way, 

i- = ~2,i ( / / J (13) 

— (-UHs°)°(-Uf 

Making use of the identities in Fig. [21 we can further use 
that X K = g R oxog A -F R oxoF A = £ R o(a;-7 R oio7 A )og A 




FIG. 2: (Color online) Identities that hold between the six 
quantities T, T, Q, G, 7, and 7 as denned in Eqs. (0-(JH]l. 
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FIG. 3: (Color online) The distribution functions x and x (left) connect incoming advanced and outgoing retarded propagators. 
The Keldysh components X K , X K , y K , and y K are shown in a diagrammatic representation of Eqs. (|14[) -(117 [ ). "R" and "A" 
refer to "retarded" and "advanced" . The particle distribution function x and hole distribution function x are coherently mixed 
due to multiple coherent Andreev scattering events with amplitudes given by the renormalized quantities Q, Q, T , and T that 
are sums of terms shown in Fig. [T] 



and similarly for the other components, which gives 

X K = g R o(x-j R oioj A )og A , (14) 

y K = g R o( X oj A - 7 R oi)o^ A , (15) 

X K = g R o(x - 7 R oa;o 7 A )og A , (16) 

y K = g R o( X oj A -j R o X )og A . (17) 

The Keldysh amplitudes X, X, y, and y are shown in a 
diagrammatic representation in Fig. [3] Note that for the 
Keldysh components we need to keep track of retarded 
and advanced coherence functions. As advanced func- 
tions propagate backward in time, their group velocity 
is reversed. Advanced propagators can be described as 
usual by the particle-antiparticle paradigm, that in the 
present case is equivalent to a particle-hole transforma- 
tion as described in Appendix [Gj In drawing diagrams 
we prefer to keep the particle picture instead of intro- 
ducing antiparticles (which would reverse the arrows and 
turn them into hole propagators with opposite energy). 

We stress that there are no diagrams with more than 
one x or x, vertex, as no retarded propagator can enter 
either of them, and no advanced propagator can emerge 
from them. As a result, the structure of the equations 
for X, X, y, and y formally corresponds to that of a 
linear response with a perturbation that switches from 
retarded to advanced (in fact, the linear response theory 
for retarded and advanced coherence functions has many 
formal similarities with the Keldysh part of the transport 
theory^ see also Appendix [E| . 



1. Alternative distribution functions 

Other definitions for distribution functions have been 
introduced in the literature. We discuss the issue of the 
various possibilities in defining distribution functions and 
their relation with each other in detail in Appendix [D] 
The distribution functions h introduced by Larkin and 
Ovchinniko v 16 ! 51 and F introduced by Shelankov 23 are 



related to the distribution functions x and x by 

x = F - 7 R oFo7 A = h + 7 R o/io7 A , 

x = F-7 R oFo 7 A = h + ^oho^. (18) 

Series expansions for the inverses can be obtained by 
iteration^ for example 

F = E"=o(7 R )"°^(7 A ) n (19) 
with (...)" = (. ..)™- 1 o(...), and 

h = E^=o(7 R 7 R )"°( a; " 7 R °£°7 A )°(7 A 7 A )"- (20) 
In equilibrium, 

h cq = F cq = tanh = -F eq = -h eq (21) 

holds. The advantages of the functions x and x are that 
the transport equations take their simplest form^ their 
numerical evaluation is easier, they simplify considerably 
time-dependent problemS ) 24 ' 54 i 61 and as we will show be- 
low, they allow for an effective handling of the boundary 
conditions. 



B. Transport equations 

The central equations that govern the transport phe- 
nomena have been derived in Ref i 24 i 32 ' 54 . The transport 
equation for the coherence functions 7 (e,p F ,i?, t) and 
7 (e,p F , R, t) are given by 

(ihv F - V + 2e)7 R ' A = [ 1 oAo 1 + £o 7 - 7 oi; - A] R - A , 
(ihv F -V - 2e) 7 R ' A = [70^07 + £07 - -foS - ii] R ' A . 

(22) 

For the distribution functions x(e,p F , R,t) and 
x(e,pp, R,t) the transport equations read 

(ihv F -V + ihd t )x - [joA + S] R ox - xo[Aoj- S] A 
= - 7 R oS K o^f A + A K oj A + 1 R oA K - S K , (23) 

(ihv F -V - ihd t )x - [ZjoA + £} R ox - xo[Aoj - £} A 
= -7 R o£ K o 7 A + zA K o 7 A + 7 R ozl K - if. (24) 
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In Appendix [E] we discuss properties of the solutions of 
these equations, and equivalent formulations in terms of 
integral equations. 

An important property of the set of equations |22|) - 
(|24[) is their invariance with respect to gauge transfor- 
mations. There are two types of gauge transformations 
that are important, and that are very different in nature. 
We discuss this issue in Appendix[Fl The first type is the 
usual gauge invariance that links the phase of the coher- 
ence functions with the electromagnetic potentials. The 
second type leaves retarded and advanced quantities in- 
variant and affects only the distribution functions x and 
x and the Keldysh part of the self energies. It leads to 
a certain freedom for the choice of the distribution func- 
tions (several choices have been mentioned above). In 
particular, when a reference system is present, distribu- 
tion functions can be defined with respect to those of the 
reference system. They are then called "anomalous",— 
and vanish in the reference system. This is particularly 
useful for situations when a system is coupled to a reser- 
voir. 



using the relations of Appendix IE 11 it follows as 
7 R ' A (P) = [ih + e^ h o5 o{e^ h + c^oSoy 1 ]^ ,(28) 
with Sq' a = [70 — 7h] R A and 

C R < A (p) = [C oe lpfh -e lp ^oC a ] R ' A (29) 
where Cq ' a is the solution of the equation 

[c Q on l - n 2 oC ] R - A = A R ' A . (30) 

For equilibrium we have E RA — —E B ' A = £ R - A , and if 
£r,a an( j [^^]r,a are di a g 0na i m S pi n space, then with 
[]?- A = -Q B - A = n B - A the relation 



7 R ' A (P) 



7 /2 + ztan(p/2)(£/Yo + A) 



fl - itan(pn)(E - j A) 
follows, in agreement with Ref. I&i. 



(31) 



3. Equilibrium solution for sub-gap energies in the presence 
of an inhomogeneous order parameter in the clean limit 



1. Homogeneous equilibrium solution 

In the case that both £ R - A = e - (£ R A - lf" K )/2 and 
[Z\Z\] R,A are diagonal in spin space, the homogeneous so- 
lutions for the coherence functions in equilibrium can be 
written as, 



7h,eq 



E ±i\/-AA-£ 2 



(25) 



where the upper (lower) sign holds for the retarded (ad- 
vanced) functions. For a singlet superconductor in the 
clean limit [Z\Z\] R,A = — \A\ 2 . In the presence of a con- 
stant superflow with momentum p s one has to replace e 
by e-v F p s . 

For the distribution function in equilibrium one obtains 

Zh,e q = (l-7 R 7 A )tanh(^). (26) 
Note that <y A = (7 R ) f (see Appendix [G|. 



2. General solution for homogeneous self energies 

For homogeneous self energies we can express the so- 
lutions 7 R,A (p) along a certain trajectory with path vari- 
able p (defined by the trajectory parameterization R = 
R + pv F ), for a given initial value 7 R,A (0) = 7q' a , in 
terms of the homogeneous solution 7 R,A that satisfies 

[ 7h o^o 7ll - Eo Jh + lh oE ~ A] R - A = 0, (27) 



where E R - 



e - E R - A , E R - 



S . Defining 



If we can neglect impurity scattering, and the system 
outside the scattering region is asymptotically homoge- 
neous with gap Ah, then for sub-gap energies |e| < \Ah\ 
we can make some more general statments about the 
properties of the coherence amplitudes. In particular, 
if we e.g. consider a pure singlet pairing state, and if the 
order parameter is of the form A ~ Ao(p)e tx ia y with spa- 
tially varying modulus Aq and spatially constant phase 
X, then, using the ansatz 7 R (p) = ij (p)e^ x+<1 '^ ■ ia y 
with real j and \P, the equilibrium equations of motion 
along any fixed trajectory read 

d 7o 
dp 

7o^ = -(l + 7o)A)Sm(^)+2e7 , (33) 
dp 

where + is a positive infinitesimal. The first equation 
is stable only in direction of increasing p. Now, for the 
initial condition far away from the scatterer, for sub-gap 
energies |e| < |Z\h| the relation j = 1 holds. Then, as 
Eq. 1321 shows, this property will be preserved along the 
trajectory regardless of the spatial variation of Ao(p). 
That means, only the phase >P of the coherence amplitude 
varies, and we have 7 R = ie l ^ x+ *wl . ia y with 



(l- 7 2 )Z\ cos(^)-0+7o (32) 



dP(p) 
dp 



-2A (p)am(#(p))+2e (34) 



and initial condition "^(po) = 0. If e = 0, then the co- 
herence amplitude stays constant along the trajectory. 
Similarly, we obtain -y R = — ie~^ x+ ' p ^ ■ 



Mjy with 



dp 



= 24)00 Bin (#(p)) + 2e - 



(35) 



Q R - A = [E - 7 h oZ\] R ' A and Q^ A = [E + Z\o 7h ] R - A , and 



For energies |e| > \A\ l \ both the modulus and phase of 7 R , 
7 R vary in space. A similar consideration can be made 
for any unitary order parameter. 
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III. SCATTERING THEORY 

We consider in the following a general quantum me- 
chanical scattering problem that is characterized by in- 
coming and outgoing Bloch wave solutions. We assume 
that the scattering region is localized in a certain space 
area, where we have in mind e.g. an interface, a sur- 
face, or an impurity. In quasiclassical context there will 
be trajectories that enter and leave the scattering region. 
Correspondingly we can define incoming solutions along 
each trajectory as those for which the group velocity is 
pointing towards the scattering region under considera- 
tion (the "scatterer" ) , and outgoing those for which the 
group velocity is pointing away. The projection of the 
group velocity on the Fermi momentum has one and the 
same sign for 7 R , 7 A , and x, and the opposite sign for 7 R , 
7 A , and x. Correspondingly, these six objects for each 
trajectory always group into three incoming and three 
outgoing ones. 

The scatterer will lead to a mixing between the tra- 
jectories that enter the scattering region. Depending on 
symmetry constraints, the possible number of scattering 
wave vectors might be drastically reduced, as for example 
is the case for conservation of parallel momentum at an 
atomically clean interface. In the latter case, for a single 
Fermi surface on each side of the interface, there will be 
mixing only between the incoming, reflected, transmit- 
ted trajectory, and a fourth trajectory that is reached by 
a process involving "crossed" Andreev reflection. In the 
case of a diffusive interface trajectories of all directions 
will be mixed with each other. 

In order to distinguish incoming and outgoing direc- 
tions we will adopt the notation of Ref. 32, that small 
case letters 7 R,A , 7 R,A , x, x denote incoming quantities, 
and capital case letters r n A , F ' , X, X denote outgoing 
quantities. As the quasiclassical Green's function is pa- 
rameterized by the momentum p F , it is composed of both 
incoming and outgoing quantities. We may write the 
Keldysh matrix Green's function as a functional of the 
four coherence functions and the two distribution func- 
tions. If the Fermi velocity points towards the scatterer, 
this functional dependence will be 

g^g[ 1 R ,r R ,r\j\x,X], (36) 

and for the case that the Fermi velocity points away from 
the scatterer, it is 

g = g[r^r,l\r A ,X,x]. (37) 

Usually the potentials in a scattering region vary on a 
energy scale large compared to the superconducting gap 
or the temperature. In this case, it is sufficient to know 
the normal state scattering matrices for particle-like ex- 
citations, denoted by S with elements S(p F <— p' F ), and 
for hole-like excitations, denoted by S with elements 
S(jj' f <— p F ), that connect outgoing with incoming quasi- 
particles on trajectories parameterized by the Fermi mo- 



menta p F and p' F J£. The scattering matrix in particle- 
hole space reads 

In order to reduce the amount of notation we will in the 
following label trajectories with the Fermi velocity point- 
ing away from the scatterer simply by k, k' , k± etc, and 
trajectories with the Fermi velocity pointing towards the 
scatterer by p, p' , p\ etc, thus omitting the vector nota- 
tion. It is understood that those labels are from the set of 
Fermi momenta associated with all the trajectories that 
overlap with the scattering region. As for the discussion 
in this chapter the dynamical variables (energy, time) en- 
ter only as parameters, we will suppress the dependence 
on these. We assume for the scattering problem that 
the spatial coordinate R on each trajectory entering or 
leaving the scattering region is sufficiently far from the 
scatterer in order that the scattered waves have taken 
their asymptotic form, but sufficiently close to neglect 
spatial variations on the scale of the coherence length in 
the scattering region, and we will suppress these spatial 
coordinates in this chapter as well. The scattering prob- 
lem will thus be fully characterized by the set of k and 
p values associated with all involved trajectories. In a 
centro-symmetric system (or a non-centrosymmetric sys- 
tem with time reversal symmetry) , for each k value there 
will also be the trajectory with the opposite direction 
p = —k. 

It is our task to express the set of outgoing coherence 
and distribution functions r^,r p , r A ,r fe , Xk X p by the 
incoming ones 7p ,7 R , 7 A ,7p , x p Xk for a given scattering 
matrix Sk P (the scattering matrix S p k for hole-like exci- 
tations is related to that for particle-like excitations by 
the particle- hole conjugation symmetry). 



A. Elementary interface Andreev scattering events 

The central objects for the formulation of boundary 
conditions for the coherence functions and distribution 
functions are the following quantities, that express an 
elementary scattering event, 

Wkk>T = E% T P °^'] H , (39) 
p 

[% P '] A = [J2 S pWk°S kp >] A , (40) 
fe 

X W = S^ p ox p oS pk ,, (41) 
v 



7 



s = 



■ ■ p ~ 
/ 3 = 



h 




x = • 



A 



FIG. 4: (Color online) Diagrammatic symbols for the ele- 
mentary scattering events described by Eqs. (|47[) - (148[) . p and 
h refers to "particle" and "hole", and "R" and "A" to "re- 
tarded" and "advanced". A sum over internal variables ac- 
cording to Eqs. (J3nj)-(I1T]) is implied. 



together with the respective particle-hole conjugated 
quantities, 



k 



pp 



(42) 
(43) 
(44) 



As we will show below, the scattering matrices enter the 
boundary conditions only in terms of these quantities. 
This allows for a compact matrix notation. For example, 
we can re-formulate boundary conditions for spin active 
interfaces that are known in literature ; 34 ' 35 in a rather 
compact way. Importantly, a straightforward general- 
ization of these boundary conditions to multiple bands, 
to disordered interfaces, to strongly spin polarized ferro- 
magnets, to strongly spin-orbit split bands, and to the 
general scattering problem from a target is possible. For 
equilibrium we recover also the results by Shelankov and 
Ozana,— that were obtained by a similar procedure. In 
order to switch to a compact matrix notation, we intro- 
duce the diagonal matrices 

llw = 7k s kk', Jpp'=JpSp P ', x kk ,=x k 5 kk >, (45) 

lpp> = lp <W > Iw = l k <W , Zpp' = x p 5 PP > . (46) 

With this, we can write the elementary scattering events 
as 



7? 



[ So 7 o £f ' A , x' = S R o xo S A 



(47) 



[ 7 ] R ' A = [So 70 S] R - A , x' = S oxo S . (48) 

In Fig. H] we show these scattering events in diagram- 
matic form. We note that the retarded and advanced 
scattering matrices are related by fundamental symme- 
try 

S A = [ S R ]t, S A = [ S R }\ (49) 

which leads together with the symmetries in Appendix [Gl 
to the symmetry relations 



x'=[x']\ 
x' = \ x'l t . 



(50) 
(51) 




fc' 



FIG. 5: (Color online) Diagrammatic representation of 
Eq. fl 53 p for the retarded functions. In the last line the iden- 
tity (|56[) is shown diagrammatically, which defines the dia- 
grammatic expansion for Q R . Summation over internal vari- 
ables is implied. 



B. Derivation of boundary conditions 

1. Retarded Propagator 

The anomalous functions T R are obtained from a sum 
over all virtual multiple Andreev scattering events that 
are accompanied by interface scattering. We consider 
first the set of retarded Green's functions with directions 
k that are directed away from the scatterer. In this case, 
the retarded propagator is given by 



g R k =mr R ,r)- 



(52) 



We introduce effective interface coherence amplitudes as 
solutions of the equation 



{Tkk-Y 



ikk> +J2f kkl oj kl o 7 ' kik , 



(53) 



Using a compact matrix notation, the solutions are 



F R = 7 'o ( 1 - 70 7 ')~ 



(54) 



where the inversion Q~ is defined via Qo Q 1 = 1 
with l kk i = 5 kk : 1 . The diagrammatic representation of 
this expansion is shown in Fig. [5l We also define the 
corresponding particle-hole diagonal interface amplitude 

Q R = [( 1- I'nrf^ + I^f, (55) 

that is closely related to the function F R by 

T R = [ go if . (56) 

For the quasiclassical approximation only the component 
k' = k is relevant, as it contributes to the slowly varying 
envelope function of trajectory k, and we obtain 



y k — y kk i •> k — ■'ki 



A- A- 



(57) 



The remaining two retarded Green's function matrix 
components are Q k = [( 1 — 70 7 / )~ 1 ]/cfc and T k = 
Gu°Jk- According to section lTl Al for the outgoing coher- 
ence functions the equation T kk = r k + T kk o^ k or k = 
(1 + ^kk°lk)°^k holds, which according to Eq. [55] is 
equal to Q kk °r k . Thus, we extract the outgoing coher- 
ence amplitudes from the solution of the equation 

Qkk° r k^k' = -Ffek'j r k = r k^k ■ (58) 

The more general quantity r k <_ k , that is introduced here 
will be needed below, e.g. in the transport equations for 
the distribution functions. 

For the component r we must consider the retarded 
Green's function for a direction p towards the scatterer, 
as the group velocity of r is opposite to the direction of 
the momentum. The corresponding retarded propagator 
is given by, 

g* = %h R ,r R ]. (59) 

We obtain in complete analogy to the discussion above 



T = 



7'o 1 - 



-A-l 
70 7 ) 



7'o 7) 



1 



Qo 7 
To 7 



(60) 
(61) 



from which we extract the outgoing coherence amplitude 
by solving the equation 

C R oT R , = F R , f R =f R (62) 
^pp p< — p pp ? p p< — p ' \ J 

2. Advanced Propagator 

For the advanced functions we need to take into ac- 
count that they propagate backward in time. Thus, we 
consider for r A the advanced Green's function for a di- 
rection p towards the scatterer, 

#=#[^,7*], (63) 

and for P for a direction k away from the scatterer, 

g£=g£h A ,r A }. (64) 

The most convenient form of the corresponding equations 
is 



T A = 7'o ( 1 - 70 7') = 7 'o g 
g A = [(1-70 7') -1 ] A = 1 + [ 7° ?\ 

with the the coherence amplitudes 



j-iA 'A <£"A nA TnA 

p' — >p ^ pp ''p'pl p p — >p ! 



and 



T = 

g A = 



■y'o ( 1 — 70 7') 



( 1 - 70 7') 1 = I 



; i'o g] 

70 T 



with the the coherence amplitudes 

^k'^k°S A k = -^fe'/s! r k = r k ^ k 



(65) 
(66) 

(67) 

(68) 
(69) 

(70) 





b) 




FIG. 6: (Color online) Diagrammatic representation of 
Eqs. J72]) and [f?5]l. using the identities in Eq. J56J and (|60)l . 



3. Keldysh Propagator 

The corresponding expressions for the Keldysh com- 
ponents are obtained in a similar way. We perform a 
diagrammatic expansion of the Keldysh components in 
the elementary scattering events, using the fact that the 
vertices x and x can only occur once in each diagram 
(see end of section Hi A[) . Thus, all renormalizations af- 
fect only the particle-hole conversion processes on either 
side of the x and x vertices. We consider first the Keldysh 
Green's function for k being directed away from the scat- 
terer, 

k K = ^ K [r R ,f\7 A ,r A ,A,i], (71) 

for which the expansion, shown in Fig. [6^, gives 

X K = g B o( x 1 - 7 ,R o xo 7 ,A )o g A . (72) 
We obtain the distribution functions X in terms of X K 

by 



Gf k o [X k - no Xk or k j oQ£ k = X^ k . (73) 

Similarly, considering the Keldysh Green's function for p 
being directed towards the scatterer, 

5 p K = 5 p K [7 R ,r R ,r A ,7 A ,x,A], (74) 
we obtain from the expansion shown in Fig. [5Jd 

X K = g B o( x' - 7 ,R o xo 7 ' A )o g A , (75) 
and from this X in terms of X , 

g; p o (x p - r;o Xp or A ) og- p = x« . (76) 
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4- Boundary conditions for coherence amplitudes 

For the outgoing coherence amplitudes that where ob- 
tained in Eqs. (53]), (1)7]), and (JO]), closed equa- 
tions can be derived, that can again be represented di- 
agrammatically in a straightforward way. We can cast 

^kt-k' = (&kk) 1 °J 7 kk'' r p ^ p , — (Gp P ) lo ^pp', rp^ p = 

FpA^p)' 1 ' and P «^k = ^ASt k r\ * the form of 
Dyson-type equations, 

[A^] R = [i kk ,+ £ ^ fel o7 fcl o 7 [,^] R (77) 

k^k 

[r P ^ P '] R = [% p , + ^p^vi°1pi % 1P 'T ( 78 ) 

and 

[^ fe ] A - [7^ fc + E 7fe^ 1 °7fe 1 °^ 1 -fc] A • (80) 

From those we obtain the quasiclassical coherence ampli- 
tudes, 

— ; = r k ^ k . (82) 

The diagrammatic representation of these equations is 
the same as for the functions T R ' A and JF RA , with the 

I 



modification that in all internal sums the direction k of 
the final state that is scattered into is excluded. 



5. Boundary conditions for distribution functions 

Analogously to the discussion for the coherence ampli- 
tudes we derive now the boundary conditions for the dis- 
tribution functions. For this we formally solve Eqs. ([73]) 
and (f76]) . 

x k - r£oz k or£ = [(SfcfcrT^^K^rT (83) 
x p - rloxpotp = [{gppy^oXppo^gpp)- 1 ^ (84) 

and use the relations 

[(£ fefe r 1 °£fcfc'] R = 1- [7V fc ,o 7fc ,] R , (85) 
[Gp>p°{Gp P r l ] h = l-l%>oT p ,^p] A , (86) 

and the corresponding particle-hole conjugated equa- 
tions. In these equations we have introduced the scat- 
tering parts of the coherence functions, that is obtained 
by subtracting the forward scattering contributions, 

[7V fc ,] R = [r k ^ k , - r k 5 kk ,] R (87) 
[r P '^ p ] A = [Tpi^p — r p Sppi] A (88) 

and similarly for particle-hole conjugated quantities. 
Solving Eqs. (|83|) and (|84|) for X k and X p leads to an 
explicit solution in terms of x' k , x k , i' p , and x p , 



X k = fa* + r *^*i 7*i] a °<fe,°[*k a fc + lk 2 or k2 ^ k ] A - ^[7V fcl ] R o£ fel o[f M ] A (89) 

fel,fe2 fcl 

x p = [ § ppi + ^P-Pi°7pi] Ro Sp 1P2 o [^2P + iP2 o r P2 ^p\ A - ^[r p^pi\ R ° x Pi°[ r Pi^p\ A ( 90 ) 

Pl,P2 Pi 



The diagrammatic representation of these equations is 
the same as for the functions X Ti A and X I< A , with 
the modification that in all internal sums over virtual 
particle-hole or hole-particle conversion processes the di- 
rection k of the state that is scattered into is excluded. 
The scattering into the final state (forward scattering) is 
taking place only in the last scattering event. Note that 
these simple diagrammatic rules result from our particu- 
lar choice of the distribution functions. Applying a gauge 
transformation of the type discussed in Appendix IF 21 to 
the distribution functions amounts to shifting terms be- 
tween the two contributions on the right hand sides in 
Fig. [6] back and forth. This leads to redefined distri- 
bution functions without changing the Keldysh Green's 
function. 



6. General use of boundary conditions 

Equations (|39 ]) -([4l ]l . ([7? ]) -([82 ]) and (JHZj) - (EOJ) give the 
outgoing quantities r£, r pl r A , r kl X k , X p in terms 
of the incoming quantities 7 R , 7 R , 7 A , 7 A , x pi x kl and 
are the main result of this paper. For a small number of 
trajectories involved in the scattering process these equa- 
tions can be solved analytically. For numerical calcula- 
tions, in particular when many trajectories are involved 
that mix with each other in the scattering region (diffu- 
sive scattering), it might be of advantage to use matrix 
algebra and solve the set of equations (|4Tl)- d48j) , (f54"]) - (f5"8]) . 
(ETJ-dSa), &CH}]), C3l-CZ3|l. and C3])- OS). The solu- 
tion of these equations involves only standard numerical 
linear algebra and is straightforward. 
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IV. APPLICATION I: SPIN-ACTIVE 
INTERFACE SCATTERING IN 
SUPERCONDUCTING DEVICES 

In this section we show how to recover from our formu- 
lation of boundary conditions the results of Refs. [H, [H], 
and [35I . These boundary conditions are for an interface 
between two superconductors or two metals or one super- 
conductor and one metal. On both sides of the interface 
each trajectory is doubly degenerate due to the spin de- 
gree of freedom. The interface is assumed to conserve the 
momentum component parallel to the interface, pn. It is 
assumed that only one Fermi surface sheet is present in 
each material, such that only one incoming and one out- 
going trajectory exists for each side of the interface. For 
such a case the boundary conditions take a particular 
simple form. As on either side of the interface (index 1 
and 2) only one incoming and one outgoing momentum 
direction are coupled by the boundary condition, we can 
label the involved trajectories simply by indices 1 and 
2, and incoming and outgoing components by small and 
capital letters in the boundary condition. 

We start with writing down Eqs. (f3"9"| - (flTj) for this case: 



7n 7i2 

721 722 



Su D12 

S2I 522 



71 
72 



521 £22 



(91) 



and 



X91 3)9 



oil S12 

521 522 



xi \ / 5n 5i2 1 , ... 

° X 2 ° 521 5 22 ) ^ 



Introducing these into the corresponding 11- and 22- 
components, i.e. Eqs. (JHSJ) and and analogously for 
the advanced functions, gives the first set of boundary 
conditions for the coherence functions, 

[A] R ' A = [7n+7i2 (l-72°722r lo 72°72i] R ' A ( 101 ) 
[AP A - [722 + 72i°(l-7i°7iir 1 °7i°7i 2 ] R ' A - (102) 

The particle-hole conjugated equations are obtained by 
simply applying the particle-hole conjugation operation 
to these results. These boundary conditions, together 
with the definitions (f9"Tl) . are equivalent to the bound- 
ary conditions of Ref. [34 ^nd for spin-scalar scattering 
matrices to those of Ref. [32|. 



B. Distribution functions 

Turning to the the Keldysh components, we formulate 
Eq. [55] for our case, 



2^1 



+r^ 2 o(i2°^2°i2 - ^)°/^l (103) 



X 2 = 4 2 + ^i°7i R| 



"12 



' 21 



>7foJV 



+AVi°(7r°^'ii°7i A - ii)°^ 2 (104) 



Substituting Eqs. (|^7T> - (|100|) into these gives the re- 
quired boundary conditions for the distribution func- 
tions. Again, the particle-hole conjugated equations are 
obtained by simply applying the particle-hole conjuga- 
tion operation to these results. These boundary condi- 
tions, together with the definitions (|92p . are equivalent 
to the ones of Ref. [H, and for spin-scalar scattering ma- 
trices to those of Ref. [33. 



where all involved quantities are 2x2 spin matrices. 



A. Coherence functions 

Using these quantities, the boundary condition 
Eq. ([77]) takes on the form 

[A^i] R = [7u + A^2°72°7 2 i] R ( 93 ) 

[A^ 2 ] R = [712 + A^2°72°7 22 ] R ( 94 ) 

[A^i] R = [721 + A^i° 7 i°7ii] R (95) 

[A^ 2 ] R = [7 2 2 + A^i°7i°7i 2 ] R - (96) 

The equations for the 12- and 21-components, Eqs. (19"4")) 
and (|95p. can be solved directly, 

[A^ 2 ] R = [7i2°(l - Wm)" 1 ]* (97) 

[A^i] R = [72i°(l-7i°7n)-T- (98) 

Analogously, for the advanced components Eq. (|79")1 leads 
to 

[A^ 2 ] A = [(l-7u°7i) _ Vi2] A (99) 

[A^i] A = [(l- 722°72)- 1 o 72 i] A - (100) 



C. Spin-active interface in bilayer geometry 

As an application we discuss the coherence functions 
for a bilayer that consists of a thick superconducting layer 
(that we will treat as bulk system) with a thin normal 
metal overlayer of thickness d. We consider a spin-active 
interface with a scattering matrix 

5= ( " s tsN )=§*. (105) 

We assume that the interface has a unique quantization 
axis, in which case all reflection (rsf, rsj., fNTi r Ni) and 
transmission amplitudes (£snt 7 tsNl, tNSh ^nsj.) are s P m 
diagonal. We consider a singlet superconductor with (re- 
tarded) coherence amplitudes 7 R = ^sicjy. As a result, 
all possible induced correlations in the normal metal are 
written as 7 R = diag[7N|, 7Njicr y , where 'diag' denotes a 
diagonal spin matrix with the diagonal elements as indi- 
cated. In the following we restrict our discussion to the 
equilibrium situation. 
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Eqs. (fT0Tj) - (fT02|) result into 



r NT r Nl7NT +*NST i SNl7S 



+7S7S7Nt(?-NT?-ST +^NST^SNT)(''Ni r Sl +*NSJ.*SNl) / 



/ 1 + 7s(^snt^ns|7n,t + r ST r s|7s) 



(106) 



Now, using the unitarity condition of the scattering 



i,>i. 



matrix, we write (with a = {T,l}), ^No- 
r S * = r.e***', t NSa = t tr e i &»*' l t s ... 
where i3 S a + i?n<t = ^sn<t + ^NSo-, and r 2 a + t\ = 1. 
Then, with the spin mixing angles $g — ~ ^SJ, and 
i?N = ^nt — ^nx, and with the further abbreviations 
[(tfsN T +tfsN;)-(tfNST+tfNS;)]/2 - (i?n+^s)/2 - 0+, 
(i?N — $s)/2 = the last equation becomes 



r- 



-ii9 4 



7n|7S' 



*T*i + 7s7se" 

i9' 



7NT e l, '-r T r i + 7se -" , 't T tj r + 7N T 7s7se lt ' + (107) 



The extra spin-scalar phase may appear due to time 
reversal symmetry breaking by the interface. In order to 
obtain the coherence amplitudes at the outer surface of 
the normal layer, 7bt, we solve the transport equation 
(ihv-p ■ V + 2e)ji(x) — in the normal metal with perfect 
reflection at the outer boundary, which gives 



7bt = r^^e 



iz/e d 



7nt = 7B| e 



iz/e 



(108) 



with z = e + i0 + , Ed = hv Fx /2d the ballistic Thou- 
less energy, and v Fx the Fermi velocity component nor- 
mal to the interface in the normal conductor. We now 
concentrate on sub-gap energies, |e| < \A\. Substitut- 
ing Eq. (fTOBf into Eq. (fTUT)) . and using the bulk solu- 
tions 7s = ie 1 ^ — — 7 with W = arcsin(e/Z\) (see sec- 
tion [TlB^J, we obtain the following equation for Jbi-> 

(109) 



7b>^' 



+ 27Bte —— 



1 



with u T = sin + + W) + r T r x sin + i?_ - 
For 7bx an analogous equation holds, with the quantity 
U[ = sin - i? + + + r T r| sin (f- - i?_ - S'). Fi- 
nally, for the particle-hole conjugated coherence ampli- 



tude one obtains 7BX e 
amplitude is given by 



-it?' 



-7bt 



AW 



Thus, the pairing 



/b<t = -27ri 



7b,, 



^T*l e " 



1 



-,2 2W 
7B.cr e; 



for \u a \ < t|t|, and by 



/bct 



iT*l si S n ( u o 



(110) 



(111) 



04 -(^i) 2 

for | Tier I > tj^x- This characteristic spin dependence of 
the pairing correlations has been discussed recently in 
Ref. 184 . where it was shown that a change in the symme- 
try of the pairing correlations near the chemical potential 
takes place as function of $n- A more detailed discussion 
will be provided in a future publication. 85 



V. APPLICATION II: 
SUPERCONDUCTOR/HALF-METAL HYBRID 
STRUCTURE 

A. Interface scattering matrix 

Next, we consider as application an interface between a 
superconductor and a completely polarized ferromagnet, 
a half metal, in the ballistic limit. Each trajectory in the 
superconductor has a spin degeneracy, whereas in the 
half metal the spin for each trajectory is fixed. Following 
Ref. l37l we write the scattering matrix in singular value 
decomposition as 



o _ ( U s 
' b ~ 1 U F 



?s \t) 
(t\ -rp 



vj 



where 1 1) and (t | are the transmission amplitudes, and 
fs = yl— \ t)(t\ and r F = y/1 — (t \ ■ \ t) are the reflec- 
tion amplitudes. The phase-matrices on the left and 
on the right can be written as Us = e l (^v+^rn u cr) ^ 



m 



a il(>u 



t _ 



singular values are determined by the matrices 



rs 
'f 



r 
1 



l*) = 
(t\=(t 



and the 

(112) 
(113) 



with r = vl — t 2 . The quantization axis is the direction 
of the magnetization in the half metal, M, which we 
chose as the z-axis. The directions are determined by 
the interface properties, and do not necessarily coincide 
with that of the half metal. 

We now make the simplifying model assumption m u = 
m v = m. We write m x = sin a cos d>. m v = sin a sin d>, 



cos a, and for the bulk magnetization M z 



M, 

M x = My = 0. Because we consider singlet superconduc- 
tors we have the freedom to choose a spin quantization 
axis inside the superconductors in a convenient way. The 
most convenient choice is along the interface magnetic 
moment m. The spin rotation matrix between the quan- 
tization axis in the superconductor and in the half metal 
is U m = e ~ i ^ e±CT with e_L = (m x M)/(M sma). In 
this representation UmUsU^ — e l ^ La * and UmVgU^ — 
e i— a z become spin-diagonal. Because in quasiclassical 
approximation only the envelope of the wave is relevant, 
we are furthermore allowed to drop all spin-independent 
phases in the scattering matrix (except for a possible 
phase ■&' analogous to that in the last subsection, arising 
from an internal flux; one can prove that all other spin- 
scalar phases do not enter the final expressions). This 
leads to the scattering matrix in the new framed 



Rs \T) 
(T\ -R F 



(114) 



e 2 







-u m 
1 



\t) 

(t\ -rp 



ut 



7 * 
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B. Josephson geometry 

The Josephson effect in a superconductor/half- 
metal/superconductor (S/HM/S) junction has 
been studied previously both experimentally 66 and 
theoretically . 37,41 - 42,55 - 67,68 - 69,70,71,72 - 73,74,75,76,77 - 78,79 - 80 
Here we demonstrate how the present formulation of 
boundary conditions can be used to simplify analytical 
expressions within the same approximation as in Ref. [76|. 
Our formulation is in terms of the microscopic scattering 
matrix. Such a scattering matrix cannot in general be 
obtained by solving Eilenberger's equations but must 
be obtained by a full microscopic quantum mechanical 
treatment of the interface ^ 8 ' 70 This has to be contrasted 
to the case considered e.g. in Ref. [7q where an interface 
represented by a thin magnetic domain wall is treated 
with Eilenberger's equations. The two approaches are 
complementary and have non-overlapping ranges of 
applicability. 



Coherence amplitudes 



solution is then given by 



r F = 



Kx-4>) 



where we use the abbreviations 







Pt 2 (l 

J2 



C = 1 



-7s7s) 
■ 27S7S?" cos($) - 
7s7s r2 + 27s7 S rcosO?) 



2-2 

7s7s 



7S7S-P 2 * 4 



(123) 

(124) 
(125) 
(126) 



assuming that all incoming coherence amplitudes at the 
superconducting side are singlet. The full solutions of 
the boundary conditions in the superconductor can also 
be obtained analytically and are given in Appendix [H] 
Note that the geometric angle </> that determines the di- 
rection of the interface magnetic moments enters only in 
combination with the superconducting order parameter 
phases. Thus, it leads to simple shifts in the current 
phase relation i 70 i 71 In the following, we include cj> into 
renormalized superconducting phases x m order to sim- 
plify notation, i.e. we define x'i — Xi ~ 4>i, X2 = X2 — §2 
for the two superconducting banks (indices 1 and 2). 



We express the boundary condition it in terms of the 
matrices 



7 " I <7l' 7f ' ^ 



7s 
7f 



(115) 



with 7s being a 2x2 spin matrix and 7f a scalar, and 
similar notations for the particle-hole conjugated com- 
ponents: 7s and 7f- Explicitely, 

% = i?s7s4+ \T) lF {f\, (116) 

1 7)' = Rsis \f)- \T) lF R F , (117) 

( 7 |' = (T\%R s -R FlF (f\, (118) 

7 F = (T| % \f)+R FlF R F . (119) 

Then the boundary conditions, Eqs. ((77)) and (|5Tj) . can 
be solved for ts and r F , leading to 



^s - % 



-1 7 5 , l7) , (7l 
1 - 7F7 F 



(120) 



= 7f + (7|' (i-7s7s) 7S 1 7)'- (121) 

This gives the outgoing amplitudes in terms of the in- 
coming ones. The particle-hole conjugated quantities are 

obtained similarly, with the definition 7g = S^jsS^ ■ 

We assume singlet superconducting order parameters 
A n = \A\e lx ia y , allowing us to write for the bulk coher- 
ence functions 7s = ^se lx ia y and 73 = 7se -JX z<x 
useful to introduce the parameter 



7se lx i*j y . It is 



P = sin (tf/2) sin(a)/(l + r) 



(122) 



with the spin mixing angle $ — + *& v , that controls the 
overall magnitude of the proximity effect. An analytic 



2. Josephson current 

The equations for the coherence amplitude in a point 
in the middle of the half metal of an S/HM/S junction, 
for positive (7+) and negative (7-) directions, can be ob- 
tained by expressing and Pp2 in terms of 7^1 and 
7F2 using the boundary conditions Eq. ll23l for each inter- 
face, and solving the transport equations in the half metal 
with the results 7+ = yTpx, 7^2 = U1+ , 7- = J/-Tf2) an d 

7fi = 2/7-, where y = e ""^^ ; and vp x is the component 
of the Fermi velocity in the half metal perpendicular to 
the interfaces. This leads for a symmetric setup to 



7+ 



7- 



y 



y 



ay '7_ — iPjse lXl 
C - if3y7se~ tx 'n~ 
a yi+ ~~ i(3jse lX2 



(127) 



(128) 



In principle the amplitudes 75 and 7 s must be obtained 
by solving self-consistently for the order parameter sup- 
pression near the interface. Here, we will however neglect 
this effect and assume that the bulk solution 



7s 



= -7s = i\A\/{e n + n n ), n n = VI^I 2 + 4 (129) 

is present all the way to the interface. This approxima- 
tion becomes exact in the limit of small t and 1?. Note 
that in this case 1 — 7 | = 2J?„/(|e„| + Q n ) is even and 



1+7| = 2e„/(|e„| + Q n ) is odd 
9+ -9- in f 1+ 7+7+ 



in 

y 



1 - 7+7^ 



Pi 



One obtains 
1 + 7_7- 
1 - 7-7- 



1 - 



(130) 
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by solving the equations 7^ + P±l± + 1 = 0. Here, 



1 

P± 



iPlsy (Ce 



2 +?'2 

rjy e ■ 



(C 2 -r;V)±2*(/37s2/) 2 sin( X ) 



(131) 



with x = Xa " Xi- Note that C - V = * 2 (1 ~ 7s7s)- 

Using this, one can show that for Matsubara frequen- 
cies p- = p*<_. Consequently, the Josephson current is 
given in terms of these quantities by, 



j = eN-pVp ■ 2nT Im 1 

e„>0 




where [i = cos(# p ), and 9 P is the impact angle (/i = 1 for 
normal impact). Here, vf and N-p are the Fermi velocity 
and the density of states at the Fermi level in the normal 
state of the half metal, respectively. 

We obtain the corresponding Josephson current for the 
case of the half metal replaced by a normal metal if we 



replace (3 = it 2 , rj 



+ isis, C = 1 + T 7s7s and add 



a spin degeneracy factor 2. 

The normal state boundary resistance of the symmetric 
S/HM/S Josephson junction with area A is given by 



1 



RkA 



j± 2 AT I t 
— = e jVpfp ( fi— — - , 

V \ 2 — t 2 / fs+ 

(dpp) . *(P||) 



(« F ei)>o 



WP||); 



t(P||) 2 



2-t(p,|)2 



t(P||) 2 



(133) 
(134) 

(135) 



with (dp-p) = d D ~ 1 pY j '(2-kK) for D dimensions, and 
P\\ = Pf ' e ||- 81 For an S/N/S junction an additional 
factor 2 has to be added on the right hand sides. 

In Fig. [TJa) we show results obtained with Eq. 11321 
Shown is the critical Josephson current multiplied with 
the normal state resistance obtained by Eq. 11331 For def- 
initeness we present results for identical isotropic Fermi 
surfaces on both sides of the interface and for the de- 
pendence of the transmission amplitude t on the impact 
angle 8 p (measured from the surface normal) appropriate 



for a (5-potential, t( 



t cos ( 



'1 



t\ sin 2 9 P , where 



to is the transmission for normal impact. For the spin 
mixing angle $ we assume a dependence 1? = $0 cos(# p ). 
For comparison we also show the corresponding val- 
ues for a superconductor/normal-metal/superconductor 
(S/N/S) Josephson junction. As can be seen, the super- 
current through a half metal can be of a similar magni- 
tude as through a normal metal, provided the parameter 
P is of order one. 

In fact, as can be seen in Fig. [TJb) , the 7 c i?N-product 
can exceed that for an analogous S/N/S junction. The 
reason for this enhancement are current carrying An- 
dreev bound states below the gap energy, that are dis- 
cussed further below. We show for several values of i9o 
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FIG. 7: (Color online) Critical Josephson current I c multi- 
plied with the normal state resistance Rn for an S/HM/S 
Josephson junction with magnetic interfaces (thick lines) and 
for an S/N/S junction with non-magnetic interfaces (thin 
lines). In (a) for both cases, the transmission amplitudes to 
are varied from 0.2 to 0.5. The spin-mixing angle for normal 
impact is $0 = 0.75-7T, and the junction length L is equal to 
the coherence length in the half metal, £f = vf/2ttT c . In (b) 
we fix L = 3£f, to = 0.5, and vary for the S/HM/S junc- 
tion the spin mixing angle. For (a) and (b) the interface spin 
misalignment angle is a = n/2. 



the / c i?N-product in comparison with that for a non- 
magnetic S/N/S Josephson junction with the same trans- 
mission probability and same length. With increasing $0 
the magnitude of the effect increases, and the maximum 
in the temperature dependence moves to lower tempera- 
tures. 

In fact, for the special case that P = 1 (i.e. t = 1, 
& = 7r, a = it/2) for all Fermi surface points, the max- 
imum becomes unobservable because it moves to zero 
temperature, as has been noted also in Ref. [76|. In 
this case, furthermore, we have (3 = ( = (1 — 7s7s)j 
i] = — 7s7s(l~ 7s7s) for the S/HM/S junction, and j3 = i, 
C = 1, T] = 7s7s for the S/N/S junction. Consequently, 
after canceling the common factor (1 — 7s7s) in Eq. 11311 
for the S/HM/S junction, it is seen that X/p± at phase 
X for the S/HM/S junction coincides with l/p± at phase 
(x + t) for the S/N/S junction. This proves that the 
J c i?N-product for P = 1 is equal to that for the corre- 
sponding S/N/S junction, and the corresponding current 
phase relations are shifted by ir. This result is in agree- 
ment with the findings in section III.D of Ref. [76] for the 
short and long junction limits, that were obtained within 
the more general Gor'kov formalism. 

We caution however, that the suppression of the sin- 
glet order parameter at the interface cannot be neglected 
for P close to 1, unless a strong Fermi surface mismatch 
is present (in which case the transmission is reduced due 
to the Fermi velocity mismatch), and self-consistent cal- 
culations must be performed as done in Ref. l4ll 
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3. Local density of states 

We now proceed to calculate the local density of states 
as function of energy. For this we need to perform an 
analytical continuation to the real energy axis. We define 
in this case y = e lzL ^ VF and 73 = -73 = -\A\/(z + 
i-y/|Z\| 2 — z 2 ) with z = e + iO + . The momentum resolved 
density of states is then given in the center of the half 
metal by, 



iV± 
N F 



! T R „ l + 7±7± 
—lmg ± = Re— 

TT 1 



7±7± 



Re 



1 - 



P±(e) : 



The local density of states is obtained as 



N{e) = (N+ + N-) 



FS, 



(136) 



(137) 



In Fig. [5] we compare the local density of states (LDOS) 
for an S/HM/S-junction and an S/N/S-junction in the 
high-transmission limit for a symmetric setup. For clar- 
ity of presentation we have vertically shifted the curves 
with respect to each other. The junctions are current bi- 
ased, and the phase difference varies in both cases from 
to tt as indicated. For the S/N/S-junction the well-known 
Andreev-Saint- James state a 20 i 82 are seen (for a review see 
Ref. |83T ) with a reduction of the LDOS at low bias except 
for the case of \ — tt, when a zero bias bound state is 
present. In contrast, for the S/HM/S-junction there is 
a low-energy band of bound states. This behavior has 
already been noted in Ref. |4l|. Note that x — 71 i s the 
equilibrium phase of the S/HM/S junctional The dis- 
persion of the Andreev peaks in the spectra with x indi- 
cates the direction of the current that is carried by them. 
For the S/N/S-junction the lowest bias peak dominates, 



that carries current in positive direction, whereas for the 
S/HM/S-junction the low-energy band is responsible for 
the low-temperature anomaly J C (T), and the next higher 
band carries most of the current, that is in negative di- 
rection, in accordance with the 7r-junction behavior. 

The half- width Wi/ 2 of the low-energy band varies 
with the interface parameters, with the impact angle, 
with the phase difference Xi with temperature, and with 
junction length. In general the width of the low-energy 
band is larger for x = than for x — 7r - hi Fig- 13 we 
show its dependence on the junction length for (a) a tt- 
junction and (b) a zero-j unction. In the short-junction 
limit the half- width for t = 1 is given by W\/ 2 (x — 7r ) — 
\A\(y/2-P 2 - P)/2 and W 1/2 ( X = 0) = \A\y/l -P 2 . 
In the limit of small P (but still t = 1) we obtain 
W 1/2 ( X = tt) -> \A\/y/2 and W 1/2 ( X = 0) -> \A\. For 
the special case P = 1 the spectra are equal to those for 
an S/N/S-junction with the junction phase shifted by 7r 
(see the corresponding discussion in the last subsection), 
i.e. the low energy band vanishes in the limit P — > 1 
for a tt junction and a zero energy bound state appears 
for a zero-junction. In general, as P varies with $ and 
thus with the impact angle, the width of the low-energy 
band in Figs.[5Ja) and [5] is a superposition for many dif- 
ferent P. The overall width of the low-energy band is 
set by the values for smallest P, the kink-features closer 
to the chemical potential correspond to the largest P for 
trajectories with normal impact. 

In Fig.[Tn]we show results for the variation of the LDOS 
with with the spin-mixing angle $0 for (a) a 7r-junction 
and (b) a zero-junction. For the zero-junction a peak 
appears at high $0, that is a signature of the zero bias 
bound state for normal impact when P = 1. For smaller 
spin-mixing angles in general the structures get smeared 
out, and a set of energy bands separated by narrower 




FIG. 8: (Color online) Local density of states in the center of 
a current biased high-transmissive symmetric Josephson junc- 
tion for (a) an S/HM/S junction and (b) an S/N/S junction. 
In both cases the phase difference over the junction is varied 
from to 7r. The remaining parameters are indicated. 
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a=7i/2 

25 ~ 9 =o.8tt 




FIG. 9: (Color online) Local density of states in the center of a 
current biased high-transmissive symmetric S/HM/S Joseph- 
son junction, for (a) a phase difference over the junction of 
tt and (b) of 0. In both cases the junction length is varied 
from the short junction limit to L = 50£f. The remaining 
parameters are indicated. 
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FIG. 10: (Color online) Local density of states in the cen- 
ter of a current biased high-transmissive symmetric S/HM/S 
Josephson junction, for (a) a phase difference over the junc- 
tion of 7r and (b) of 0. In both cases the spin mixing angle t?o is 
varied from to 7r. The remaining parameters are indicated. 



suppressions of the LDOS remains. 

C. Point contact geometry 

1. Distribution functions 

For the distribution functions, we introduce the nota- 
tion, 



\ x F J 
with S A — (S R y, which explicitely gives 

x' s = R%x s R A + x F \T) R (T\ A 

\x)' = R%x s \T) A -x F R A |T) R 

(x\' = (T\ R xg R A x F R F (T | A 

x' F = (T\ R xs \ T) A + x F Rf R$. 



with Ri 



\ry< 
(rr 

and|r> A = «f |*)t, (r 
ary conditions, Eq. 
read, 



(138) 



(139) 
(140) 
(141) 
(142) 

T)*)t, 

(143) 
(144) 



= (|r) R ) t , the explicit bound- 
for the distribution functions 



an i T) A = «T | R )t, (T y 
Then, with the abbreviations 



I t) / (1 — 7f7f) _1 



(7 1' (1 -7s7s 



*s - x' s +^\xY(r\^+^\r) R (x\' 

+(^x' F ^-x F )\mr\ A (145) 

x F = x' v + (x\' % \r) A + (r\ 



a 7 S R \x) 



+ UT (7s R is7s A -^s) \h A (146) 



2. Point contact spectra 

Superconductor/half- metal point contact spectra 
have been studied experimentally in a number of 
caseSl 2&£LSS£2£&SLS2£&2i2£ However, the analysis in 
all these studies did not include the effect of spin active 
interface scattering. Here, we show how such effects 
can be taken into account in a ballistic point contact. 
We assume incoming solutions to be in equilibrium. 
The treatment in terms of coherence and distribution 
functions can be simplified considerably by using the 
symmetries described in Appendix IF 21 Proceeding along 
the lines described there, we introduce anomalous distri- 
bution functions by g K [x, x] — g K [xo, xo] — g T 'Fo ~ Fog A , 
with 



Fn = 



F 








(147) 



and Xo = x — (F + 7 R i^o7 A ). We use for F the equilib- 
rium distribution function in the superconductor. Then, 
the incoming anomalous distribution functions xs.o and 
xs,o m the superconductor are zero. For the half metal 
we have x F = F + -f F F% with F = tanh[(e - eV)/2T] 
and F = — tanh[(e + eV)/2T], and consequently xp.o 

I — F ) + 7p(F — Fq)^. Furthermore, for a ballistic 
point contact the incoming coherence functions on the 
half-metallic side are zero, 7f = 7f = 0. From here on 
we drop the index "0" for all distribution functions in or- 
der to not overload the notation, and keep in mind that 
they are all anomalous. 

Substituting all this into Eq. (|146[) . one arrives at 



X F — x F — 



x F t 2 



IC R f 

r cos(a) sm(d)A 2 



{(1 - n 2 ) [l + r 2 n 2 - r cos(i))S 2 



-t 2 P 2 {l + r) [E 2 {1 + rU 2 ) + 2(1 + r)77 2 ] } 



-r^x F f F = -^P 2 (l + r) 2 | 7s R | 2 [1 + S 2 + n 2 ] 



(148) 



with 



n 2 



l7s7sl 



(7s B 7s R )(7s A 7s A ) 2 = 
-(7s R 7s R + 7s A 7s A : 

- 7 (7s7s R " 7s A 7s A ) = -2Im( 7s R 7 s R ) , (151) 



(149) 

2Re(7s R 7s) (150) 



and we have used the notation 7 g = ^icjy and 7 g — 
7gi(jj,, meaning that 7 g = — (7s)*- These expressions 
are still general and we only made use of zero incoming 
7 f, 7f, x s ,o, and x s ,o- 

The current density from the half metal to the super- 
conductor is given in terms of the anomalous distribution 
functions by, 



Here,7 s A = (7 s R )t, 7 - = (7f) 



j = eN F v F 



% L{X F x F r$x F f F )\ (152) 
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FIG. 11: (Color online) Point contact spectra for an S/HM 
contact. In (a) the transmission to varied from to 1, and in 
(b) the spin-mixing angle $o is varied from to tt. Both quan- 
tities depend on the quasiparticle impact angle as discussed 
in the text. The remaining parameters are indicated. 



For comparison, we also present the expressions for the 
corresponding di/d^-spectra for a normal metal, 

„2|i /-„,n\2|2l 



|1-- 2 (7 S R )T) FS+ 



Je = 

This 



and 



l + (Mi 4 l7sT-r 2 |i-(7s R ) 2 

gives j e = <2 M < 4 /|l-r 2 ( 7s «) 2 | 2 ) Fg+ , for |e| < \A\, 

U = <M* 2 [1 + (7s R ) 2 ]/[l - ^ 2 (7s R ) 2 ]) FS+ for |e| > \A\. 

In Fig. [TT]we show representative results for zero tem- 
perature S/HM point-contact spectra for various trans- 
missions and spin- mixing angles. In general, there are 
sub-gap states present in the spectra except for very small 
■do, a, or t . For the special case do = it there is a sharp 
zero bias state observable in the spectra. Otherwise, if 
i9o is not close to tt, dl/dV is for T = zero at zero bias 
and increases quadratic with the voltage. The details of 
the spectra will depend on the Fermi surface mismatch 
and the interface characteristics, in particular the depen- 
dence of the various parameters on impact angle. We 
leave a detailed discussion of these issues for a future 
publication. 



with x-p = F — Fq = [tanh e — tanh ^?], and x F = 
F — Fo = — [tanh e \^f — tanh Wy]. The current density 



can be written as 

'°° de 



j = eN F v F 



oo 



tanh ■ 



eV 



2T 



tanh ■ 



eV 



2T 



(153) 



with spectral current kernels j e . The normal state bound- 
ary resistance is 1/R^A = jw/V = e 2 iV F v F (fJ-t 2 ) fs+- 

Further simplifications arise for incoming homogeneous 
distribution functions, when 7 g ~ — 7 g. Then, noting 
that for |e| < we have Hi = 1, it follows that 



Xp — xp 



-r F i F r£ 



Xpt 
X-pt 



P 2 (l + r) 2 {E 2 + 2) (154) 



P 2 (l + r) 2 (S 2 + 2) (155) 



leading to the Andreev spectral current 

/ ^-4P 2 t 4 (l + r) 2 • (l + Re[( 7 g) 2 ]) 



\ 1 + (7 R ) 4 r 2 - 2( 7 <?) 2 r cos(tf) + ( 7 f ) 2 P 2 t 4 I 

FS + 

(156) 

For zero misalignment of the interface moments, a = 
leads to P = 0, and there is no Andreev current. For 
|e| > \A\, additional terms become important, associated 
with n 2 1. Here, because 7 g is real, and thus i F (e) = 
— a; F (— e), we obtain 

M* 2 ■ [1 ~ (7s R ) 4 

1 + (7s R ) 4 r 2 - 2( 7 £) 2 rcos(tf) + ( 7s R ) 2 P 2 *\ FS 

M .2P 2 t 4 (l + r) 2 .[l + ( 7 f) 2 ] 2 ( 7s -) 2 \ 



1 + (7s ) 4 ^ 2 - 2( 7g , ) 2 rcos(tf) + ( 7 £) 2 P 2 t 4 



FS 4 



(157) 
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APPENDIX A: TIME CONVOLUTION 
PRODUCT 

We use extensively the non-commutative o-product be- 
tween two functions, which allows us to formulate the 
equations independently from the representation of the 
dynamical coordinates (time, energy, mixed). In the 
time domain, the non-commutative o-product between 
two functions A(t,t') and B(t,t') is defined by 



AoB(t,t') = / dt" A(t,t")B(t",t'), 



(Al) 



with the unit element 1 = 5{t — t')l. In an energy repre- 
sentation (after a Fourier transform t — > e, t' — > e'), the 
product reads 

AoB(e, e') = J ^ A(e, e")B(e", e'), (A2) 

with the unit element 1 = S(e— e')l. In a mixed represen- 
tation, when performing a Fourier transform (t~t f ) — > e, 
and keeping the time variable (t + 1') jl — > t, the product 
can be written as 

AoB(e,t)=e^ 9 ^- 9 ^A(e,t)B(e,t) , (A3) 



17 



and the unit element is i = 1 . If one of the factors is both 
independent of e and t, the o-product reduces to the usual 
matrix product. Note that in a mixed representation 



eoa — aoe = ihd t a, eoa + aoe = 2ea. 



(A4) 



Sometimes (for example when performing a perturbation 
theory out of the equilibrium) a modified energy repre- 
sentation is useful, where one performs Fourier trans- 
forms (t — t') — > e, (t + t')/2 — > oj. In this case the 
product reads 



Analogously, for the linear response to an exter- 
nal perturbation, the normalization condition leads to 

-R.A pR.A = 5; as a con- 



P 



~Sg R ' A oP*' A = and 
sequence the spectral response, <S? R,A , can be written as 



Sg R - A = ^2iti P+oSWoP_ - P_o5ZoP + 



(B4) 



with a suitable parameterization of the functions SW and 
SZ. 



AoB(e,w) 



du>' du>'' 



6(u)' + uj" -u)x 



2ir 2?r 

xA(e+^,iv")B(e-^f,co') , (A5) 



and the unit element is 1 = S(u>)l. If A(e,t) = A(e) is 
independent of t (if A is an equilibrium quantity) then 



AoB(e,u) = A(e+!f)B(e,u) , 
and, analogously, if B is an equilibrium quantity 



(A6) 



AoB(e, u) = A(e, w) j?( e _ *2) . (A7) 
We also generalize throughout the paper the commutator 
[A,B} = AoB-BoA . (A8) 
A useful identity is 

(I + aoby 1 oa = ao (1 + boay 1 . (A9) 

APPENDIX B: PROJECTORS 

We adopt here the notation of Ref. H2, Appendix B. 
Following Shclankov^ we introduce the following pro- 
jectors 



(Bl) 



with the properties P + oP + = P + , P_oP_ = P_, P + + 
P- = 1, and P+oP_ = P_oP+ = 0. The quasiclassical 
Green's function is expressed in terms of P+ or P_ by, 



= -in (P + - P 



(B2) 



From the normalization condition, the Keldysh compo- 
nent of the Green's function, g K , fulfills the relations 
P^og K oP A = and P R og K oP A = 0, which allows a pa- 
rameterization by 



9 



-2?r i 



P%X K oP A + P*oY K oP A 



(B3) 



were X K and Y K are related by symmetry relations. The 
function X K can be chosen in a convenient way. 



APPENDIX C: PARAMETER 
REPRESENTATIONS OF PROJECTORS 

The projectors Pf and P R may be parameterized by 
complex spin matrices 7 R and 7" as defined in Appendix 
C of Ref. l32j. Alternatively, we give here the parameteri- 
zation in terms of Q, J 7 , G, and T. We obtain 



PI 



Pi = 



9 T 

-t (l-g) 



-T 
9 



,P* 



(l-g) -t 

T G 



G T 
-f (1-9) 



(Cl) 



(C2) 



APPENDIX D: PARAMETER 
REPRESENTATIONS OF DISTRIBUTION 
FUNCTIONS 



In general the functions X K and Y K in Eq. (|B3[) can 
be written as 



X K 



Xll 


%12 \ 




J/11 


J/12 \ 


2/12 


Vn J 




I X12 


in J 



(Dl) 



taking into account the fundamental symmetry relations 
for the Keldysh Green's function through the "tilde" 
particle-hole symmetry relation. Any choice of the four 
functions xu, X12, j/n, and j/12, will lead to a valid pa- 
rameterization of the Keldysh Green's function. As they 
parameterize only one free function in g K (due to sym- 
metry relations and normalization condition), three of 
the four parameters can be chosen conveniently. It is 
customary to require X12 — j/12 = 0, leading to the pa- 
rameterization 



X K 



x 
y 



y 
x 



(D2) 



Three definitions for distribution functions have been 
considered in literature. They correspond to different 
choices of the remaining two parameters. Larkin and 
Ovchinnikov introduced the parameterizatio n 16 ! 51 



h : X K 



h 
-h 



,Y K 



-h 
h 



.(D3) 
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Shelankov's distribution functions 23 follow from 

F 



x = y = F : X K = 



F 
F 



F 

The author introduced the parameterization 24 
y = y = 0:X K = ( I ° V Y K 





x 



(D4) 



(D5) 



The advantage of (|D5|) is that the transport equations 
take their simplest form. The advantage of (|D4[) is that 
X K and y K are scalar in particle-hole space. And the 
advantage of ([D3]) is that X K + Y K = 0. Why the latter 
property is an advantage one can see when re-writing 
Eq. (B3J) into 



.9 = 



ITT 

y 



(X K + Y K ) + — o(l K - f K ) - (X K - f K )o — — o(X K + Y~ K )o — 

— in —in —in —in 



(D6) 



r 



With X K + Y K = this leads to 
g K = g R oX K - X K og A with X K = 



h 
-ft 



(D7) 



which is an equivalent definition to Eq. (|D3[) that was 
first given by Larkin and Ovchinnikov. 

The symmetry relations for all these distribution func- 
tions are 



h(e,p F ,R,t) = h(-e,-p F ,R,t)*, 
F(e,p F ,R,t) = F(-e,-p F ,R,t)* 
x(e,p F ,R,t) = x(-e,-p F ,R,t)*, 



and 



h(e,p F ,R,t) 
F(e,p F ,R, t) 
x{e,p F ,R, t) 



h(e,p F ,R,tf, 
F(e,p F ,R,tf, 
x(e,p F ,R,tf. 



(D8) 
(D9) 
(D10) 



(Dll) 
(D12) 
(D13) 



The x and x are expressed in terms of the other distri- 
bution functions in a straightforward way, and we obtain 
Eqs. (fIS ]) -([2D j) of the main text. 

Finally we comment on the linear response, Eq. (|B4[) . 
Here, the most convenient parameterization is 



Sy R - 




5Z B 







(D14) 



APPENDIX E: PROPERTIES OF THE 
EQUATIONS OF MOTION 

In this Appendix, we use some shorthand notation in 
order not to be confused by too cumbersome expressions. 
We use for the superscripts (R, A, M) the notation (X). 
We parameterize the position on the trajectory by a spa- 
tial coordinate R = pv F . We also introduce the symbol 
d for hv F ■ V and omit the o symbol in all products. 
Finally, we use the shorthand notation E x = e — S x , 
E x = -e- £ X , and E K = -S K , E K = -if. 



1. Relations between different solutions for 
coherence functions 

We consider solutions of the equations of motion for 
the coherence functions, Eq. 1221 and for simplicity we 
concentrate on the first one, as the second is related to the 
first by fundamental symmetry relations. The equation, 



■ J 



A 



-r + £ X 7 X - 7 X ^ X + ^ x = o, 

7 X (0) = 7 f, (El) 

is a Riccati matrix differential equation, the basic prop- 
erties of which were thoroughly studied, e.g. in the book 
of ReidiS 2 - Associated with any solution 7 x (p) of (|E1[) 
are three quantities U x (p\j x ), V x (p\j x ), and W x (p|7 x ), 
which obey the set of equations 



i8U x + {E x - 7 X Z\ )U X = 0, 
tdV x -V X {E X + A x ~f x ) = 0, 
idW x + V X A X U X = 0, 



C/ x (0) = 1 , 
^ x (0) = 1 , 
W x {0) = 



(E2) 
(E3) 
(E4) 



Let us assume we know the solution 7q (/?) with initial 
condition 7* (0) = 7^ and associated functions Uq, VqS 
and Wq. Often it is the case that we have boundary 
conditions, which have to be fulfilled for given molecu- 
lar fields, external fields, and order parameters. Then we 
have to find the initial value 7^ self consistently. A prop- 
erty of Riccati differential equations is that the knowledge 
of one solution allows to construct any other solution. For 
this we note that the solutions U x (p), V x (p), and W x (p) 
for any other initial condition j x = 7^ + Sf are 

U X (P) = U x (p)[l + 8 x W x (p)}-\ 

V x (p) = [l + W ( f(p)S x ]- 1 V x ( P ), 

W x (p) = [l + W^p)S x }- 1 W^(p) 1 

= W^m+dfW^p)}- 1 . (E5) 

The full solution j x (p) along the entire trajectory for the 
new initial condition is then obtained by the following 
formula: 



7 X (P) 



7 X (P) 
7 X (P) 



UZ(p)S x V x (p) 
U x {p)5 x V,f{p). 



(E6) 
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2. Integral equation for coherence amplitudes 

For the retarded and advanced coherence amplitudes 
there is a possibility to formulate the Riccati differential 
equation as an integral equation. The formal solutions of 
equations (|E2HE4|) arc 

U*(p) = 7V ; /o(£ x -^ x )V 

V x (p) = p e -*/ r f(£ x +^ X 7 x )dp" (E7) 

V^py 1 = 7V ; /o P (£ x +^VW' (E8) 



where V ('P)is a trajectory (anti-) path-ordering opera- 
tor. With the definition of the transfer operators 

SX(p,p>) = U*( P )U*( P r 1 = Ve l ^^ A * )Ap " 



sX( P ',p) = v*(p')- 1 v x (p) = ^/;'(^ X +^ X 7 X )V 

(E9) 

and introducing the notation 

I A (p) = -A*(p) - ^( P )A x (p)^(p) (E10) 
we can write the equation of motion as 

id 7 X + (£ - 7 ^) X 7 X - 7 X (£ + ^7) x - I a (EH) 
and obtain an integral equation for 7 X , 

7 x (p) = ^(p,0)7 x (0)^(0,p) 

- i [ P S$(p,p')lZ(p')S$(p',p) dp' .(E12) 
Jo 

3. Construction of solutions for distribution 
functions 

In a similar way we can obtain integral representations 
for the Keldysh Green's functions. Consider the trans- 
port equation for the distribution function x, 

idx + (E- -/A) R x - x(E + A-/) A = I K 

I K = 1 R E K ^y A + A K ^y A + j R A K + E K . (E13) 

The solutions can be written in terms of Sfj{p,p') and 
S A {p,p'), Eq. jE9}, as 

x(p) = S^(p,0)x(0)S A (0,p) 

- i f S^(p,p r )I K (p / )S^(p',p)dp / (E14) 
Jo 



4. Construction of solutions for linear response 
functions 



Analogously we obtain the linear response equations 
for retarded and advanced coherence functions, that are 



given by the solutions of 

id <5 7 x + (E- 7 ^) x <5 7 x - Sj x (E + Aj) x = SI X 
« x = 7 X (5^ X 7 X - <5£ X 7 X + -f x SE x - SA X . (E15) 

Its solutions can be written in terms of (p, p') and 
S$(p,p'), Eq. (TE9]), as 

^(p) = S^{p,0)Sf (0)^(0, p) 

- i [" SUP,P')SI X (P')SUP',PW ■ (E16) 
Jo 



APPENDIX F: GENERALIZED GAUGE 
TRANSFORMATIONS 

We start with the set of quasiclassical equations 

[e-h,g] +idg = gog = -^\. (Fl) 

We note that the generalized gauge transformation in 
combined Keldysh and Nambu-Gor'kov space 



g' = f- 1 ogof 



(F2) 



leaves equations (|F1[) invariant, 

[e-h',g']o + idg' = , g' o g> = -tt 2 1 . (F3) 
if we use as gauge transformed source terms 

e-h! = f- 1 o(e-h)of + f- 1 oidf . (F4) 
Here, the matrix T is of the following form 



T = 
We write now 

f = f D +f K = 



T A 



/ fu 


V 






V 






o J 



(F5) 



(F6) 



Here, the matrix To is assumed to have an inverse T D l . 
Then, the inverse of f is expressed through Tp by 
f- 1 = T D X - T D X o f K o f D \ Defining f K = -f D o F 
we can write without loss of generality 

f = f D o(l-F) f- 1 = (l + F)of D 1 (F7) 

where the matrix structure of F is given by 

F 



F = 







(F8) 



and F o F = ensures the simple structure of the inverse 
of T. Now we can write the generalized gauge transfor- 
mation as 

g' = (l + F)of D 1 ogof D o(l- F) (F9) 
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Thus, we have two types of transformation, which we can 
study separately, first for F = 0, 

g' = fjj 1 ogof D 
e — h! = t D x o (I -h)of D + f' 1 o idf D . (F10) 

and second for To = 1, 

g' = (l + F)ogo(i-F)=g-[g,F} , 
e-h' = {e-h)-[e-h,F] -idF . (Fll) 

The second transformation does not affect the re- 
tarded and advanced components, and redefines only the 
Keldysh components. It leads to a gauge transformation 
for the distribution functions. A general transformation 
is obtained by successive application of these two types 
of transformations. 

For an infinitesimal transformation Td = 1 — STd we 
obtain to first order 

g' = (1 + 5f D ) o g o (I - ST D ) = g- [g, 5T D } , 
e-h' = {e-h)-[e-h,6f D ] -idSf D . (F12) 

Note the similarity to the second type of gauge transfor- 
mation. This follows from the fact that formally we can 
write (i + F)= e p , and (1-F) = e~ p due to FoF = 0, 
so the same equations like for Td = e~ 5Tr > w 1 — STd 
hold. 

1. Transformations of coherence functions 

The Riccati differential equations are invariant under 
the following transformation with transformation matri- 
ces f x and T x , 



7o X = (T x )- 1 o 7 x oT x , (F13) 

A x = (f^^o^or, (F14) 

Eq = (f x )- 1 o{idf x + E x of x ), (F15) 

Uo = (f x )- 1 o[/ x of x , (F16) 

Vo = (T x )- 1 oV rX oT x , (F17) 

W x = {T x )- 1 oW x of x , (F18) 

x = (f R )- 1 oxof A , (F19) 

A* = (T n )- 1 oA K oT A , (F20) 

Eq = (T R )- 1 oE K of A , (F2I) 



with E x = e — S x , E K = and analogous relations 

for the particle-hole conjugated quantities. An important 
case is that of unitary transformation matrices, where 
this transformation is a local gauge transformation, pos- 
sibly accompanied by a local spin rotation. In this case 
it is more convenient to write 

T x = f x = e~5<A . (F22) 



The important feature is the occurrence of the new driv- 
ing terms (T x ) _1 o idT x that gives a contribution to 
the vector potential. For gauge transformations they are 
equal to 

-^=e5*o^o F .V^oe-'f (F23) 

When vp-V</> commutes with <p, the two gauge factors on 
either side cancel in equilibrium. As can be seen above 
there is a very broad class of transformations (not neces- 
sarily gauge transformations) which leave the equations 
of motion invariant. 

2. Transformations of distribution functions 

The equations of motion are also invariant under the 
transformations 

x = i-(F +7 R ofo°7 A ) (F24) 
A$ = A K + (A R oF Q + F oA A ) (F25) 
E$ = E K - (E R o F - F o E A ) - idF (F26) 

with E x = —E x , E K = —S K , and analogously for xo, 
A Q , and Eq. A natural choice is the equilibrium distribu- 
tion function, F = tanh(e/2T) (and F related by sym- 
metry) . The transformed quantities are called anomalous 
in this case. 

Let us assume we calculate the Keldysh Green's function 
from x and x and obtain <? K [x,:r]. Applying the above 
transformation of the driving terms we could also solve 
for the xo and xq instead. We can then construct an 
anomalous Green's function defined by g a = g K [xo,x ]. 
The difference between the Keldysh part and the anoma- 
lous part of the Green's function is called spectral part 
of the Green's function. If we introduce 

«■ = ( -ft ) ™ 

then it is given by, 

g K [x, x] - g K [x ,x } = g R o F - F g A (F28) 

Thus, it is enough to solve for xq and xq to obtain di- 
rectly the full Keldysh Green's functions once one has the 
retarded and advanced ones. The choice of the distribu- 
tion function _Fo is of course somewhat arbitrary, but it 
is best chosen to be the equilibrium distribution function 
whenever there is one well defined. For a spatially vary- 
ing electrochemical potential and possibly varying 
temperature, 

where $(R) is determined by the unit trace of the 
Keldysh Green's function to ensure local charge neutral- 
ity. The advantage of such a choice is that the anomalous 



21 



functions Xq are zero in 'reservoir' regions. If the elec- 
trochemical potentials are different on the two sides of 
an interface, then the boundary conditions produces a 
nonzero anomalous component xq on either side of the 
interface. It is always numerically advisable to use the 
xq with the spectral part subtracted instead using the 
full x. This makes the driving forces explicit and avoids 
cancellations between large terms. 

Let us finally mention the driving terms for the above 
choice of equilibrium function. They are given in the 
equation for xq by —idF , with 



dF = v 



eE(R) - Vfi{R) - 



e$(R) 



VT(R) 



hd e F 



T(R) 

(F29) 

where E is the electric field. This corresponds to the 
force term in a Boltzmann equation. There are additional 
terms for time dependent forces. For instance the term 
e o F — Fq o e is equal to ihdtFo. Note also the term 
— A R o Fo — Fq o A a which gives for energy independent 
gap as off-diagonal force 



A ■ ( tanh [ - , 



e$(R) \ , / < 

x ' 1 - tanh 



e$(R) 



2k B T(R) 



(F30) 



Finally, we mention the possibility to define spin depen- 
dent forces in a similar way. 



APPENDIX G: RETARDED- ADVANCED 
SYMMETRIES AND KELDYSH SYMMETRIES 

The following symmetries connect retarded and ad- 
vanced functions and express symmetries in the Keldysh 
components: 

7 A = (7 R ) f , A A = -(Ay, E A = (E R )i, (Gl) 
U A = (V R )\V A = (U R )\W A = (W R )\ (G2) 
x = (x)\ A K = (A K )\ P K = -(F K ) t (G3) 

with E RA = e - r R ' A , E K = -S K . The quantities U R - A , 
V R ' A , and W R - A are defined in Eqs. (|E5j) . Analogous 
relations hold for the particle-hole conjugated quantities. 



APPENDIX H: FULL SOLUTIONS IN 
SUPERCONDUCTOR FOR S/HM INTERFACE 

The solutions of the boundary conditions for the 
model discussed in section IV B II can be obtained also 
in the superconductor explicitely. With the abbrevia- 
tions ?9fj = ($ u — #„)/4, P = sinfsina/(l + r) and 
Q = cos f sina/(l + r), they are given by 



1 - t Pt 2 (l + r)[7 F7 | e ^+7 F e^]+7s(l-7 F 7F)[2rcosz9~P 2 t 4 ] 
Us)(u-lT)/2 2 l-r 2 7 F7F + iPi 2 (l + r)7 F 7 S e-^ [ ' 

1 -Qt 2 (l + r) [ 7F7 2 e -^ +7Fe ^]+ j 7s (l-7 F7F ) [2r sin + PQi 4 ] 

y 1 Sj(Tl+lT)/2 — o 1 2 : — i r m27i i — =i5 

2 1 — H 7f7f + iPi z (l + rj7 F 7ge 1<P 

, r s _ 2itf TT ~* 2 [^h'^ sin2 f ~ 7 F e^ cos 2 f ] + iPt 2 ls [(1 + r 7F7F ) sin 2 f + (r + 7f7f ) cos 2 f ] 

\* SJtT e e 2 •r>j.2/'i i —7(4 I, 

1 — r z 7 F 7 F + iFt z (l + r)-fFjse 9 

, r x - 2 « tfTT ^ [7F7s 2 e-^ cos 2 f - Tpe** sin 2 f ] + zP£ 2 7s [(1 + r 7F7F ) cos 2 f + (r + 7f7f ) sin 2 f ] 

Us)u = e rT e r - — : — ^ToTi i — F =7^ ■ ^ H4 J 

1 — r^7 F 7 F + iPt z (l + r)7 F 7se %<p 



For zero misalignment angle a, these solutions sim- (is)TT = t 2 7F e 24,,TT /(l — r 2 7F7F ), and (is) 44 = 
plify: (P S )( U -iT)/2 = 7 srcos(i?)(l- 7F7F )/(l-r 2 7F7F ), -t 2 j Fl % e -™n /(l - r 2 7F7F ). 
( r s)(Tl+iT)/2 = i 7 s?-sin(?9)(l - 7f7f)/(1 - r 2 7 F 7 F ), 
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